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The NOAA Geostationary Operational Environmental Satellites (GOES) fire 
thrusters to remain within a 1° longitude-latitude box and to dump accumulated 
angular momentum. In the past, maneuvers have disrupted GOES imaging due 
to attitude transients and the loss of orbit knowledge. If the R-series of 
spacecraft to be launched starting in 2012 were to follow current practice, 
maneuvers would still fail to meet Image Navigation and Registration (INR) 
specifications during and after thruster firings. Although maneuvers and 
recovery take only one percent of spacecraft lifetime, they sometimes come at 
inopportune times, such as hurricane season, when coverage is critical. 

To alleviate this problem, thruster firings small enough not to affect imaging are 
being considered. Eliminating post-maneuver recovery periods increases 
availability and facilitates autonomous operation. Frequent maneuvers also 
reduce longitude/latitude variation and allow satellite co-location. Improved 
orbit observations come from a high-altitude GPS receiver, and improved 
attitude control comes from thruster torque compensation. This paper reviews 
the effects of thruster firings on position knowledge and pointing control and 
suggests that low-thrust bums plus GPS and feedforward control offer a less 
disruptive approach to GOES-R stationkeeping and momentum management. 


INTRODUCTION 

The GOES spacecraft fire thrusters to maintain station and dump angular momentum. On 
GOES N-P, daily momentum dumps and bi-monthly east-west (longitude) maneuvers 
will be done during a daily 10-minute housekeeping period. Large north-south 
(inclination) maneuvers are done once a year, and INR performance may take 6 hours to 
recover. By performing small maneuvers on a daily basis, one could certainly reduce box 
size but might also reduce availability because of the time to recover position knowledge 
and stabilize attitude. This paper attempts to assess the feasibility of daily maneuvers in 
the light of expected orbit estimation and attitude control performance. 8 
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BACKGROUND 


GOES spacecraft are stationed at 75° and 135° west longitude. Without maneuvers, both 
would oscillate about the -105° longitude due to the Earth’s nonspherical gravitational 
field. Their orbital inclinations would also increase by about 1° per year. To keep the 
spacecraft close to their nominal longitudes in the equatorial plane, thruster firings are 
necessary. Tangential thruster firings lower or raise the semi-major axis to set the 
spacecraft drifting east or west. Out-of-plane thruster firings cancel inclination drift 
caused by solar and lunar gravity. In addition, thrusters are needed to dump angular 
momentum due to solar radiation pressure. 

As shown in Table 1, north-south maneuvers require about 50 m/s velocity change (Av) 
each year and east- west maneuvers require about 1 .3 m/s Av each year. If maneuvers 
were done every day, the daily Av’s would be about 13 cm/s and 0.35 cm/s. Momentum 
dumps of about 25 Nms are also needed each day for spacecraft with unbalanced solar 
arrays. It is difficult for GOES to have balanced solar arrays because of the need for a 
clear field-of-view to space for the instrument cryogenic radiators. For thrusters with 1 m 
moment arms, the momentum dump Av is about 1 cm/s. Daily north-south bums are 
largest followed by momentum dumps and finally east-west bums. 


Table 1 GOES Maneuver Schedule 



thruster 

axis 

magnitude 
Av (m/s) 

period 

(days) 

daily 
Av (m/s) 

east-west stationkeeping 

tangential 

0.27-0.29 

-80 

0.0035 

momentum dump 

any 

0.01 

1 

0.01 

north-south station keeping 

out-of-plane 

45-50 | 

365 | 

0. 13 


Synchronous Elements 

Because of the ambiguity as to right ascension of the ascending node Q and argument of 
perigee to , Keplerian elements are not well-suited for work with geosynchronous orbits. 
A commonly used alternative is the set of synchronous orbital elements. 1 Semi-major 
axis a is replaced with drift D 


D = -3Sa/2A (1) 

where A is the nominal geostationary semi-major axis and 8 a is the deviation from it. 
Eccentricity becomes a two-dimensional vector with magnitude equal to the Keplerian 
eccentricity but pointing toward perigee 
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( 2 ) 


e T = {e x e y ) = e (cos(Q + co) sin(Q + to)) 

Inclination becomes a two-dimensional vector with magnitude equal to the Keplerian 
inclination pointing along the orbit normal projection in the equatorial plane 

i T = ( i x iy ) = i (sin D - cos Q) (3) 

The eccentricity and inclination vectors are shown in Figure 1. Mean anomaly is replaced 
by sidereal angle measured from the geocentric equatorial x-axis 

s — G + X — G 0 + y/ • (t — 1 0 ) + X (4) 

where G is the Greenwich sidereal (or hour) angle, X is the sub-satellite longitude and 
\|/ is the Earth’s sidereal rotation rate, i.e. 360.985647°/day. 



Figure 1 Orbital Elements 


Position and Velocity 

From the synchronous elements, the east-west (A,), north-south (0 ) and radial (r ) 
coordinates of the spacecraft may be computed as 

0 =-*' cos 5 — i sin 5 

x y 

r = A- A(3D/2 + e x cos s + e y sins) (5) 

A = A 0 +D(s -s 0 ) + 2e x sins-2e v cos s 


Here, the subscript o denotes epoch position, i.e. where the spacecraft started out. As can 
be seen from the first equation, latitude variation is equal to the inclination magnitude. 
Similarly, the last equation shows that longitude variation is equal to twice the 
eccentricity magnitude. 
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Thruster bums may be used to control drift, eccentricity and inclination according to the 
following equations 9 

Av r = v(Ae x sinsj -Ae y coss^) 

Av, == v(aD + 2 Ae x cos s + 2Ae y sins b ) (6) 

A v 0 = v (A i x sin s b -A l cos s b ) 

where the subscripts r , t and o indicate the radial, tangential (increasing longitude) and 
orthogonal (increasing latitude) directions. The subscript b indicates that the sidereal 
angle is at the thruster burn time. 

Geosynchronous stationkeeping strategies are described in Soop and in journal articles. 
Among them are the articles by Kamel on maximizing the time between longitude 
maneuvers and on eccentricity control and a more recent one by Kelly on combining 
longitude and eccentricity control. 2 ' 4 In these papers, maneuvers are treated as 
instantaneous because commonly used chemical thrusters are so powerful that bums do 
not take much time. More recently, papers have addressed low-thrust stationkeeping. 5 ' 6 


ENVIRONMENTAL EFFECTS 

Even in geosynchronous orbits, spacecraft may move with respect to the ground. This is 
because of perturbations from the Earth’s nonspherical gravitational field, solar radiation 
pressure and lunisolar gravity. The nonspherical geopotential causes spacecraft longitude 
to drift. Solar radiation pressure drives orbit eccentricity and applies a torque on the 
spacecraft. Lunar and solar gravity cause orbital inclination to drift. These effects are all 
countered by thruster firings. 


Nonspherical Geopotential 

GOES-East and -West are stationed on either side of a geopotential minimum at 105.3° 
west longitude. Without stationkeeping maneuvers, the spacecraft longitude would 
oscillate about this point with a period of more than 2 years. 7 Thrusters are fired when 
the spacecraft reaches the “lower” edge such that the spacecraft travels “up” across the 
box and turns around at the other edge. This is depicted in Figure 2 for a 0.1° box, and 
the maneuver period is about 25 days. Longitudinal acceleration over the box is about 
0.001 “/day 2 which makes the box width for daily maneuvers 0.00025° or 4 prad. This 
corresponds to 140 meters east-west movement of the sub-satellite point. 
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Figure 2 Longitude Stationkeeping 


Solar Radiation Pressure - Eccentricity 

Eccentricity is driven by solar radiation pressure. As solar radiation pressure accelerates 
the spacecraft on one side of the orbit, the radius stretches out on the other side. There, 
the solar radiation pressure slows the spacecraft making the orbit still more elongated. 
Integrated over a day, the average rate of change in the eccentricity vector is 


de __ 3 p s A f- sins s ' 
dt 2vm y cos s s y 


( 7 ) 


where p s is the solar radiation pressure, A is the spacecraft area, m is the spacecraft 
mass, v is the spacecraft velocity and s s is the sidereal angle of the Sun in the equatorial 
plane as shown in Figure 3. 



Figure 3 Solar Radiation Pressure and Eccentricity 
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For ballistic coefficients (mass-to-area ratio) over 50 kg/m 2 , eccentricity grows less than 
0.4x1 O' 6 per day and can be controlled by choosing the maneuver time. Because 
longitude variation is twice the eccentricity, this eccentricity causes less than 1 farad 
longitude variation. GOES I-M eccentricity is currently controlled by choice of maneuver 
time, and the typical eccentricity is less than 0.0005 which gives a longitude variation of 
only 17.5 prad. Figure 4 shows the almost circular path of the eccentricity vector over the 
course of a year. If this level of eccentricity variation is acceptable, the trajectory can be 
centered on the origin with relatively small bums. 


0.0005 



Figure 4 Annual Eccentricity Drift 


Solar Radiation Pressure - Angular Momentum 

Spacecraft angular momentum increases due to solar radiation pressure and offsets 
between the spacecraft center of pressure and center of mass. For an axially-symmetric 
spacecraft with a single solar array, the primary torque component is in the orbit plane 
perpendicular to the Sun line. To compute the momentum build-up, one can integrate 
Euler’s equation in coordinates rotating with the spacecraft at angular velocity co 

dH 

M= + axH (8) 

dt 

Here M represents the solar radiation torque and H represents the spacecraft angular 
momentum. H is actually the sum of the spacecraft rigid body angular momentum fin 
and the wheel angular momentum H w 


H = Ia> + H w 


(9) 
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For a perfect attitude control system with no angular acceleration, the differential 
equation for the wheel angular momentum is 

H w =M-3x{l3 + il w ) (10) 

Because it is tied to the Sun line, the angular momentum rotates in spacecraft body 
coordinates. To estimate the size of the momentum buildup, we used GOES I-M 
dimensions minus the solar sail which balances the torque from the solar array. 8 
Assuming that the spacecraft center-of-mass is at the center of the spacecraft body, the 
angular momentum builds up in inertial space as shown at the top of Figure 5. In 
spacecraft body coordinates, it cycles between the in-plane spacecraft axes, with growing 
amplitude. Over 24 hours, the angular momentum increases by about 20 Nms. 



Figure 5 Angular Momentum Growth 


Lunisolar Gravity 


Orbit inclination is the deviation of the orbit normal from to the celestial north pole and is 
driven by solar and lunar gravity. Because the Earth is also accelerated by the Sun and 
Moon, the spacecraft acceleration relative to the Earth is written as 






( 11 ) 


where F and r k are the spacecraft and Sun/Moon position vectors in Earth-centered 
coordinates, and p and /u k are the Earth and Sun/Moon gravitational constants. 9 
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Because the inclination vector is parallel to the angular momentum vector, its direction 
obeys Euler’s equation 

^ = ? x 4t O 2 ) 

dX 

where h is the spacecraft angular momentum taken about the Earth-spacecraft center of 
mass and normalized by the spacecraft mass. In terms of h , the components of the 
inclination vector are 


i = 


' sm' (hjhf 
v sin - l (h y /h\ 


(13) 


Both Sun and Moon push the inclination vector toward the vernal equinox, and although 
the Sun’s pull on the spacecraft is larger than that of the Moon, the Moon actually has the 
greater influence because it is the differential force that drives inclination. Figure 6 
shows the annual 1° inclination drift. 



Figure 6 Annual Inclination Drift 


MANEUVER RECOVERY 

There are two reasons for the INR maneuver recovery period. First, because of Av 
uncertainty, the new orbit is not well known until there is enough post-bum tracking data. 
GOES-R 3a position uncertainty must be less than 100 m on each axis. Second, when 
thrusters are fired, the solar array may vibrate and cause attitude oscillations. GOES-R 
3 a attitude and angular velocity errors must be less than 360 prad and 100 prad/s on each 
axis respectively. Because the instruments are only required to meet their INR 
requirements once the spacecraft does, GOES-R must recover quickly. 
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Orbit Knowledge 


GPS has long been seen as the key to onboard orbit determination and rapid maneuver 
recovery. 10,11 Until receivers sensitive enough to track GPS sidelobe signals were 
developed, geosynchronous coverage was spotty, and the time needed to converge was 
long. With new receivers having 22-25 dB-Hz sensitivity, it is possible to track multiple 
GPS vehicles from geostationary altitude most of the time as shown in Figure 7. ’ 



High-fidelity simulations have been done showing that the GOES-R position knowledge 
requirement of 100 m per axis can be met during normal operations. 14 This study 
assumed the GPS Block HA antenna pattern and 9 dB receiving antenna gain. To 
estimate the time required for orbit knowledge to converge following a maneuver, a 
simple covariance analysis was performed in which the initial velocity covariance was 
increased from its steady-state value by a value corresponding to the Av uncertainty. 

The covariance P was propagated using the orbit state transition matrix O and the 
process noise covariance Q 

P<-OPO r +Q (14) 

The information matrix P~' was updated by adding the information in each observation 

P- 1 <-P~' +G t WG (15) 

where W was the observation weighting matrix and G was the derivative of the GPS 
range observations with respect to spacecraft position. 
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Figure 8 shows the position knowledge uncertainty (3a) following north-south maneuvers 
with different Av uncertainties. The observation interval is 10 seconds, the observation 
standard deviation is 1 m, and there are 4 GPS vehicles in view. With daily maneuvers, 
the expected Av uncertainty would be less than 1 cm/s. Such small uncertainties are 
practically undetectable. The solid curve is for a 1 m/s Av uncertainty comparable to that 
of an annual 50 m/s north-south maneuver. The dash-dot middle curve is for a 10 cm/s 
Av uncertainty. For GOES-R, where position is required only to 100 m, it should be 
possible to fly through daily maneuvers without violating orbit knowledge requirements. 



Figure 8 Post-Maneuver Position Knowledge 


Attitude Control 

Thruster firings affect attitude control rather than knowledge. To quantify their effect, the 
spacecraft was modeled a? shown in Figure 9. The main body was replaced with a sphere 
free to rotate through the angle (p about a diameter perpendicular to the page. Its 2000 kg 
mass and 1.5 m radius were chosen to be close to those of GOES I-M. The solar array 
was replaced by a point mass at the end of a cantilever beam. Its 80 kg mass, 5 m length, 
0.25 Hz natural frequency o) n and 0.01 damping ratio q were also chosen to match the 
GOES I-M solar array. 



Figure 9 Spacecraft-Solar Array Model 
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For small angular displacements c|> and small linear displacements 5 , the open-loop 
equations for the spacecraft and the solar array are 


(m l +m 2 )x l =F x (17) 

(m l + m 2 )y { + m 2 (R + l)<j> + m 2 d = F y (18) 

J 0 + m 2 (R + £) • ( i>, + 5) = M r + M c (19) 

m 2 ( jq + (R + £)<j> + 5) = —cS - kS (20) 


where F x and F y are thruster forces, M T and M c are thruster and control torques and 
J is a combined moment of inertia 


/ = /, +m 2 (R + £) 2 (21) 

The solar array damping and spring constants c and k are derived from the damping ratio 
and natural frequency as 

c = 2qco n m 2 — 0.5924 kg/ s (22) 

k = co 2 n m 2 = 46.5281 kg Is 2 (23) 

A proportional-integral-derivative (PID) control system as shown in Figure 10 was 
chosen because an integral term was needed to eliminate the steady-state error to the long 
thruster pulse. The control system roots were chosen to lie in a Butterworth pattern with 
natural frequency equal to the 0.01 Hz bandwidth co c . 

K r =2o C J 

K p = 2 a> 2 J (24) 

K^colJ 

The bandwidth was chosen to be well below the solar array vibration frequencies so as 
not to excite them. Knowledge of the applied torque was assumed accurate to 5% and 
was used in feedforward fashion to reduce the input torque M . 


Thrusters 


J Force 


[Torque — 



Spacecraft Attitude 


Figure 10 Attitude Control System 
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To impart 20 Nms of angular momentum to the spacecraft with a 0.1 N thrust and a 1 m 
moment arm, one would fire for 200 seconds. Figure 11 shows the attitude error and 
angular velocity predicted for such a firing. The attitude error peaks at the start and end 
of the pulse but stays below the 360 prad requirement. The angular velocity also stays 
well below the 100 prad/s requirement. If controller bandwidth could be doubled to 0.02 
Hz, the peak errors would drop by a factor of 4 to 50 prad. This is still well below the 
solar array natural frequency. Although not shown here, the response to a pure force, 
rather than a pure torque, was also computed but was much smaller. Thus, it seems quite 
plausible that attitude control can be maintained through thruster firings. 
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Figure 11 Thruster-Induced Attitude/Angular Velocity Error 


CONCLUSION 

This paper has examined the feasibility of doing GOES-R stationkeeping and momentum 
dumping with frequent low-thrust bums assuming a high-sensitivity GPS receiver and 
attitude control system with thruster torque compensation. The results suggest that: 

1. Longitude and latitude control can be tightened from the current 1° to 0.1° or less. 

2. Orbit knowledge requirements can be met during and after daily thruster firings. 

3. Attitude control requirements can be met during and after daily thruster firing. 

While encouraging, these results are preliminary, and concerns remain as to choice of 
thruster, development of a flight-qualified high-sensitivity GPS receiver and design of the 
feedforward attitude control system. 
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